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Abstract 

We report on the status of primordial nucleosynthesis in light of recent results 
on CMB anisotropies from WMAP experiment. Theoretical estimates for nuclei 
abundances, along with the corresponding uncertainties, are evaluated using a new 
numerical code, where all nuclear rates usually considered have been updated using 
the most recent available data. Moreover, additional processes, neglected in previous 
calculations, have been included. The combined analysis of CMB and primordial 
nucleosynthesis prediction for Deuterium gives an effective number of relativistic 
degrees of freedom in good agreement with the simplest scenario of three non de- 
generate neutrinos. Our findings seem to point out possible systematics affecting 
4 £fe mass fraction measurements, or the effect of exotic physics, like a slightly de- 
generate relic neutrino background. 
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1 Introduction 



The recent results of WMAP collaboration on Cosmic Microwave Background (CMB) 
anisotropics pQ perhaps represent both a change of perspective and a further call for 
increasing precision in cosmology. Before their first data release, in fact, the determina- 
tion of baryonic matter density Uf, = Qf,h 2 was in the realm of Big Bang Nucleosynthesis 
(BBN), providing a result which typically is a compromise between the values preferred 
by primordial Deuterium and 4 He experimental determinations (see for example j2j). The 
very accurate determination which is now available from CMB, u>b = 0.024 ± 0.001 [H], 
is an independent piece of information which may be used as a prior in BBN theoret- 
ical analysis. This has dramatic consequences in determining the internal consistency 
of the standard picture of light nuclei formation and/or the role of systematics in the 
experimental measurements. As we will discuss in the following, the Deuterium (density 
fraction with respect to H) experimental result of J!., X D = (2.78lg;3g) ' 10~ 5 , is in quite 
a fair agreement with the theoretical prediction using the CMB value for cj&. This is a 
remarkable success of standard BBN, since is a rapidly varying function of Uj,. The 
status of the other two observed nuclei, 4 He and 7 Li, is less satisfactory. The source for 
the tension between theory and experiments, though not so serious at the moment, may 
be due to systematics in the reconstruction of the primordial value from observations or 
rather to some more exotic features of the BBN theoretical scenario, as extra relativistic 
degrees of freedom or neutrino chemical potential. 

To clarify all these issues it is therefore necessary at this stage to make an effort to 
further improve the accuracy of theoretical predictions, as well as to have new observation 
and measurement campaigns. In the recent few years some steps have been made in this 
direction: the accuracy of the neutron to proton ratio at decoupling and the numer- 

ical stability of the original public BBN code ^0] have been improved, the uncertainties of 
the theoretical results due to experimental errors on the several nuclear rates entering in 
the reaction network ^T] have been quantified. This network, partially updated ^2] with 
the recent NACRE compilation > traces back to the original analysis of [2j and, more 
recently, of fSlEj- In the last decade, however, many new experimental and more precise 
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data have been obtained for relevant nuclear reactions in the energy range of interest for 
BBN. It is therefore timely to perform a complete and critical review of the whole nuclear 
network, in order to refine the input values for the rates and, perhaps more crucial, to 
reduce the theoretical uncertainties on nuclide yields (mainly D and 7 Li). In Section 2 we 
summarize the main results of such an analysis. All nuclear rates already present in the 
code of Wagoner-Kawano [TU] have been updated using the most recent available data. 
Moreover additional processes, neglected in previous calculations, have been included in 
a new version of our BBN code. Relevant examples of updated reactions and new added 
ones are discussed in this paper, though a complete overview is in progress [T7j . 

In the following we will be mainly interested in determining cj& and the effective number 
of relativistic degrees of freedom, N e ff, defined as 



with p 7 the photon contribution to the total relativistic energy density, Pr. In the stan- 
dard scenario three active massless, non degenerate, neutrinos are the only additional 
relativistic degrees of freedom. This results in N e jf = 3.04 at the CMB epoch, which 
takes also into account the small entropy release to neutrinos during the e + — e~ anni- 
hilation phase and QED plasma corrections PHj-[2II|- This effect also slightly modifies 
BBN and results in an effective N e ff = 3.01 for 4 He mass fraction fZT\ |2*2"1 ITS] . As known 
from many analysis |2S]-|2S1, if one considers N e ff as a fitting parameter, CMB alone 
cannot severely constrain it, though WMAP data resulted in a sensible improvement of 
previous bound. We have performed a CMB likelihood analysis, described in Section 3, 
using the CMBFAST public code [29 j and the WMAP software facilities to calculate the 
likelihood functions [SHI- [S3]- Assuming that N e ff is not changing between the BBN and 
CMB epochs, we show in Section 4 that a joint analysis of CMB data and Deuterium 
abundance nicely fits in the simplest standard scenario, with N e ff ~ 3. 

As we mentioned, the present determination of 4 He mass fraction would require a 
slightly smaller value for u>b than what is suggested by WMAP data. A possible way 
out to this discrepancy is of course to reduce the neutron to proton ratio at decoupling, 
e.g. to lower the relativistic energy density. It is therefore not a surprise that a purely 




(1.1) 
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BBN likelihood analysis of the iV e // - Ub parameters, described in Section 5, provides 
a preferred value for N e ff smaller than three. We conservatively interpret this finding 
as a possible sign of systematics in the 4 He measurements. In Sections 4 and 5 the 
estimates of primordial abundances are discussed. Deuterium results for different quasars 
absorption systems (QAS) distribute around the theoretical expected value, while there is 
still an evidence at 2-3cr level for primordial 7 Li depletion. We also report the theoretical 
prediction for 6 Li and 3 He abundances, which may represent further tests for BBN. 

In Section 6 we perform an updated analysis of the degenerate BBN scenario. Similar 
studies are reported in EH]- It has been shown in [HE] that for large mixing angle 
oscillation solution to solar neutrino problem, chemical potentials of the three active neu- 
trino species are much more severely bound than what was found in |2S] , since oscillations 
tend to homogenize the lepton asymmetries among the active neutrino flavors. In par- 
ticular, the conservative bound on the asymmetry parameters, £ a = fi Ua /T u , quoted in 
is \^fj,, T \ ~ |£e| < 0.07. Of course there is still room for a small neutrino-antineutrino 
asymmetry in the universe, which may help in reconciling the slight discrepancy in the 
A He abundance. Remarkably, a tiny neutrino chemical potential has a great impact on 
the BBN favorite value for N e fj. We finally give our conclusions in Section 7. 

2 A new code for BBN 

In this Section we briefly describe the main aspects of a new BBN code realized starting 
from the original public version of Kawano jlOj . This has been achieved in several steps, 
some being already described in previous papers [HI E]. The main improvement is a 
complete review and update of all nuclear reaction rates which enter in the light nuclide 
network, as well as, more generally, a critical review of many other processes which were 
not considered in the standard BBN network. Some of these resulted to be not negligible 
and thus have been included for the first time. An exhaustive description of the analysis 
of all rates included in the code is in progress [TZj. Here we only stress the key role of 
some of the processes we have either updated using new experimental results or included 
for the first time. 

The main aspects of the new code can be summarized as follows: 
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1) Numerical features 

As discussed in jH], numerical resolution of differential equation describing the time 
evolution of nuclide abundance is quite critical. This is because the problem is stiff, 
since the right hand sides of corresponding Boltzmann equations, being the difference 
of almost equal forward and reverse process rates, are typically a small difference of 
two large numbers. To deal with this problem our code adopts a NAG routine [3*o] 
implementing a method belonging to the class of Backward Differentiation Formulae, 
instead of a traditional Runge-Kutta solver. 

2) n <-> p weak rates and A He mass fraction 

The weak rates converting neutrons and protons are the key inputs deciding for the 
n/p ratio at freeze-out and eventually for the 4 He mass fraction. The theoretical estimate 
for these rates was improved by considering QED radiative and finite temperature effects 
as well as finite nucleon mass (Sj-jO], which result in a reduced uncertainty in A He mass 
fraction Y p at the level of 0.1%. Present statistical (and possibly systematics) uncertainty 
of the Y p experimental determination, quoted at the level of 1%, does not presently suggest 
any further effort in theoretical analysis of this issue. We have also updated the value of 
the neutron lifetime, r n = (885.7 ± 0.8)s jSZ|- 

3) An updated nuclear reaction network 

The nuclear reactions involved in BBN form quite a complicated network of order hun- 
dred different processes. The original compilation (standard network) was reviewed 
by Smith, Kawano and Malaney ^H] • The standard chain was partially reanalyzed in ^2] 
(for a more recent analysis see Ref. [TS]) by extensively using the NACRE nuclear rate 
collection. Since then new data became available, as for example the LUNA collaboration 
results for D + p — > 7 + 3 He (see later) jHH]- It is well known that the uncertainties on 
several nuclear rates still represent the main limitation to any theoretical analysis and 
prediction of BBN. On the other hand, an improved precision of BBN theoretical esti- 
mates is mandatory, both to check the consistency of the physical grounds of BBN and to 
constraint the values of cosmological parameters. This can be achieved by a new complete 
study of all reaction rate network. Our analysis can be summarized as follows: 

i) all two body strong and electromagnetic reactions entering in the production/destru- 
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ction of light nuclei up to A = 7 have been all reevaluated, using NACRE estimates 
for rates and errors. In particular the key reaction n + p 7 + D has been studied 
in details, both in its experimental and theoretical aspects, since it represents the 
main D production process, therefore influencing all the heavier nuclide primordial 
production. 

ii) for all processes not included in NACRE catalogue we have either used other avail- 
able experimental data or simple theoretical models to get a reasonable estimate for 
the rates. 

iii) new processes are now included in the code, as for example the 7 Li + 3 if <-> n + 9 Be 
stripping/pickup reaction j^H], and the analogous 7 Li + 3 He <-> p + 9 Be, 7 Be + 3 H <-> 
p+ 9 Be, which may be relevant in opening new channels for the production of heavier 
nuclei with A > 7. The rate for 3 He + 3 H <-> 7 + 6 Li has been also considered and 
added to the code. This rate was usually neglected since it is assumed that 6 Li 
is mainly produced via 4 He + D — > 7 + 6 Li, whose rate is however poorly known 
(different low energy non resonant rate estimates disagree for even one or two order 
of magnitude). 

iv) we have analyzed the possible role of some three body processes, as p+n+n <->• n+D, 
n + p + p <-> p + D, p + n + n <-> / ~f + 3 H, and n + p + p <-> 7 + 3 He. In agreement with 
standard approach, we checked that these processes can be in fact safely neglected, 
being suppressed by the high entropy per baryon during BBN. 

The benefit of this analysis is of course a higher accuracy in several rates. As a gen- 
eral statement we can say that uncertainties for 4 He and D primordial yields are now 
largely dominated by experimental errors. In principle, if no systematics were present in 
experimental determinations, a careful analysis of the consistency of the standard BBN 
scenario would be possible, providing stringent constraints on cosmological parameters uj, 
and N e ff. Further data, mainly on D — D strong reactions and 4 He + 3 He <-> 7 + 7 Be, 
are nevertheless still highly desirable. 

As an example of our study we here report some of the main new results on rates 
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Figure 1: The data and fit for the astrophysical factor of D + p — > 7 + 3 He from NACRE 
catalogue (left plot) and with LUNA data added (the fit of S(E) is with a cubic polyno- 
mial). 

included in the new BBN code. This is by no means a complete overview, which we will 
present elsewhere [T7], but it is enough to illustrate the main aspects of our analysis. 

a) n + p <-> 7 + D 

This reaction is the main process for D production. There are no many experimental 
results for the corresponding rate in the energy range relevant for BBN {3U . In 
the previous versions of the BBN numerical code the corresponding rate estimate 
[T5] was based on the measurement of the thermal capture cross-section, <j t h = 
(0.3326 ± 0.0006) barn |JT], with an overall uncertainty on the rate of 7%. 

We have evaluated the rate on the basis of the few low energy data [12] and mainly 
using the Effective Field Theory. For the low energy range, it is possible to de- 
scribe the process in terms of an expansion of local operators up to N 4 LO, with 
electromagnetic coupling introduced via gauge principle (13j. The resulting fit of 
the rate has been then smoothly linked up with the estimate of [15] . which is more 
accurate for T > 0.1 MeV. The overall uncertainty is found to be 1.2% and 7% for 
T < 0.1 MeV and T > 0.1 MeV, respectively, which represents a sensible improve- 
ment in the relevant region for BBN. 

b) D + p <-> 7 + 3 He 

This reaction is the main direct source of 3 He together with the strong process 
D+D «-> n+ 3 He. The corresponding rate evaluated in is based on experimental 
data which are accurate in the energy range 0.1 -v- 10 MeV, while for lower energies 



there are two different estimates |HJ 113 • A fit of the corresponding uncertainty 
shows a strong temperature dependence, increasing from 5% at T ~ 10 MeV up to 
30% for T ~ 0.01 MeV. 

A sensible improvement of accuracy is provided by a reduced incompatibility of 
the two set of measurements |44| l4"5j . due to the presence of systematics in [H] 
(see [IE]). More recently, the LUNA collaboration has performed a new detailed 
measurement of this process rstte £tt very low energies [38] j close to the solar Gamow 
peak, E ~ 6.5 • 10~ 3 MeV . We show in Fig. ^ the astrophysical factor as presented 
in the NACRE catalogue and the fit we have obtained including the new LUNA 
data. The rate uncertainty is now reduced to at most 3.6% (statistical error). 

c) D + D <-> n + 3 He and D + D <-> p + 3 H 

These two reactions are the dominant contribution to, respectively, 3 He and 3 H 
production for T < 0.1 MeV. Moreover they represent a main source of uncertainty 
in the theoretical determination of the deuterium final abundance (see i.e. fl]). We 
have reanalyzed the data collected in by fitting the (non resonant) astrophysical 
factor S(E) with a cubic polynomial (in good agreement with the NACRE estimate), 
but the adopted rate has been fitted with a better accuracy. In our analysis we adopt 
the NACRE uncertainty. However, an accurate evaluation of the error budget of 
the data sets, presently in progress ^7j, seems to indicate that a purely statistical 
error (below 2%) can be attained. In any case, even such a large reduction of the 
rate uncertainties has no relevant consequences on the likelihood analysis, due to 
the dominance of the experimental errors. 

d) 3 He + 3 H ^ 7 + 6 Li 

This process, modelled with an electric dipole transition [37], is usually neglected 
in the BBN network, since present estimates suggest that it contributes only for 
at most 1% to 6 Li production. This nuclide is in fact mainly synthesized via the 
electric quadrupole transition A He + D <-> 7 + 6 Li. Nevertheless, we have decided 
to include the 3 He + 3 H <-> 7 + 6 Li rate for two reasons. First of all the rate for the 
process 4 He + D <-> / y + 6 Li is only known for E > 1 MeV and around the resonance 
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at E = 0.711 MeV, while at lower energies there are only weak upper limits. On the 
other hand, the theoretical estimates still differ for order of magnitudes; thus 
its relative weight in 6 Li production may be smaller than expected. Moreover the 
3 He + 3 H «-> 7 + 6 Li rate itself is possibly affected by a normalization uncertainty, 
of the order of a factor 3 (at most) (see Ref. 8 in [I3|). Our data fit of the S factor 
for 3 He + 3 H 7 + 6 Li with a cubic polynomial shows a good agreement with the 
result of [48 . 

We define, as usual, A He mass fraction, Y p , and relative density with respect to hydrogen, 
Xt (i = D 3 He, 6 Li, 7 Li), as 

y -V A H4He 
-^AHe = ip — 4 , 

n b 

X< = — , (2.1) 

with rif, the baryon density. We do not use 7 Li in our analysis to constrain the values of 
cosmological parameters (see Section 5 for a discussion of this nuclide abundance). To 
compare the experimental values with theoretical estimates we construct, as customary 
in literature, a likelihood function £ as follows. If Q represents the set of cosmological 
parameters entering the theoretical model (basically Ub and N e ff and, in the degenerate 
BBN scenario, the asymmetry parameter for electron neutrinos, £ e ), and R the set of 
all nuclear rates Rk involved in the network, the joint x 2 {ty function (i.e. —2 logC) is 
defined by 

x 2 (Q) = E R) - XT) Fa [xf{n, r) - xf p ) , (2.2) 

ij 

where F is the inverse of the error matrix and the sum is over the nuclei considered in the 
analysis. With Xj h (Q,R) and X^ xp we denote the theoretical and experimental estimate 
for the z-nuclide abundance, respectively. In the following we will be interested in studying 
either the D or the joint D + 4 He likelihood function, 

-21og(£ D ) = (x%(n,R)-X™ p ) 2 F 22 , (2.3) 
-2\og(£ BBN ) = J2 ( X t h (^ R ) ~ X T P ) Fa [Xf- (Q, R) - Xf p ) , (2.4) 

i,j=D,AHe 

For £ D 

F£ = {<? + « P ) 2 , (2-5) 
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while for the joint likelihood D + A He the error matrix is given by 

^7 1 = (0 2 + K? P ) 2 > (2-6) 

i.e. the sum in quadrature of the (diagonal) experimental error matrix (cr^) 2 = 5ij{al xp ) 2 
and the theoretical one (cr*^) 2 (f2, R). The latter, accounting for the effect on the nuclide 
abundances of the uncertainties of the several nuclear rates R k , can be estimated using 
the linear propagation method introduced in [11]. In a completely general approach, we 
consider the quantities 

«) 2 (fi,i?) = ^[XifaRk + S^-XifaRk-SRi) 
4 k 

x [Xj(n, R k + 5R+) - Xj(n, R k - SR*)] , (2.7) 

with 5R k the (temperature dependent) upper and lower uncertainties on R k , respectively. 
This differs from the approach in ^T] which assumes the existence (in principle not neces- 
sary) of the linear functionals \i k = d\ogXi(Q)/dlogR k . This only holds for symmetric 
and temperature independent relative uncertainties 5R k /R k . 

3 WMAP results and u b 

Since we are interested in using the CMB data in a joint analysis with the BBN 
theoretical/experimental ones, the first step in this study has been to produce our own 
database and likelihood analysis of WMAP measurements. We have used the CMBFAST 
code [29, and the public software provided by the WMAP collaboration to calculate the 
likelihood function [30J- 33J. The temperature and polarization anisotropies have been 
evaluated for a set of ~ 2 ■ 10 6 models, varying the following parameters in the reported 
ranges with the corresponding steps: 
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Figure 2: The 68 and 95% CL. likelihood contours from WMAP data in the N e ff — out, 
plane. 

We only consider spatially flat models, as suggested by the position of the first acoustic 
peak, now determined by WMAP with a rather high accuracy, l p = 220.1 ± 0.8 [49] . We 
have also let as a free parameter the overall normalization, which is eventually fixed by 
the likelihood analysis. In Table Q we report the best fit values and 95% CL. ranges 
obtained after marginalization over remaining parameters. Notice that h does not appear 
in this Table since we consider its range as a prior. 

By using as fitting parameters u>b, u) c , h, r and n s and fixing N e ff = 3.04 one gets at 
the best fit x 2 = 1431 for 1342 d.o.f. . The introduction of N e ff as a new free parameter 
lowers the value of x 2 to 1430 for 1341 d.o.f., not significantly improving the goodness of 
fit. This is in perfect agreement with the results of [24 . 

The results are shown in Fig. I2HH1 After marginalization we get Ub = 0.0231qooi an d 
N eff = 2.6i?;°, at 68% C.L., in good agreement with similar results obtained in the recent 
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UJ C 
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0-981^ 


1 1 +u - uy 
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0.115 


0.023 


0.15 
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Table 1: The best fit values and the (2a) bounds for marginalized distributions are here 
reported for the fit analysis of WMAP data. 
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Figure 3: The marginalized likelihoods versus N e ff and w&. 

literature [23]-|2n]- Compare, for example, our Fig. 3 for N e ff with the red curve in Fig. 
1 of Ref . [21] . We note that the range obtained in j2H] is instead slightly larger. 

For completeness we report in Fig. 0] bidimensional likelihood contours in all planes 
involving at least one of the parameters uo^ , N e ff. From this figure we see how CMB 
data alone cannot break the degeneracies N e ff-h and N e ff-u) c , which are the main sources 
of uncertainty in the determination of N e ff. Removing these degeneracies requires other 
data sets than CMB. Notice that the ripples in the contours are an artefact of the finite 
size in the parameter lattice. 

4 CMB+BBN(Deuterium) analysis 

Deuterium number fraction, Xq, is critically depending on the baryon content of 
the universe. It is therefore worth starting our analysis of the consistency of the CMB 
+ BBN standard scenario by discussing how WMAP determination of fits with the 
experimental determination of Xr>. 

The estimate of Xp is affected by a theoretical uncertainty which is dominated by a 
small number of nuclear rate uncertainties, mainly the processes which produces either 
3 if or 3 He out of Deuterium. Varying all the rates Rk which enter the network with la 
errors 5R^ (see Section 2), we get from our code, for uo b = 0.023, cr^ — 0.15 • 10~ 5 . In 
Table |2] we also show the relative contribution of the main processes to this result, in 
percentage. 

Our previous estimate for the same parameter was er|| = 0.21 • 10~ 5 jO]. The main 
source of this improvement is due both to the more precise data which are now available 
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to Ub, cu c , h, t and n s , respectively. 
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Figure 5: The pulls (see text) of QSO D measurements with respect to the theoretical 
prediction for N eff = 3.01 and uj b = 0.023, in units of ((o^) 2 + (a e 2 xp ) 2 ) 1/2 . 

on D + p <-> 7 + 3 -£?e and the reanalysis of the D — D rates. Considering cr^ and the 
allowed range for uo b obtained in Section 3, 0.022 < u>b < 0.025, we get, for the standard 
value N e ff = 3.01, X D = (2.55±o; 2 |) ■ 10~ 5 . This theoretical prediction nicely fits, within 
the errors, with the most recent experimental results of [I], X^> = (2.78lo;3g) ' 10~ 5 . The 
latter is the average of several measurements of DI /HI column ratio in different QSO 
absorption systems at high red-shift. 

In Fig.EJwe report the deviations of the experimental values of X^, with respect to the 
theoretical prediction of standard BBN, in unit of the combined error due to cr^ and a^ v 
summed in quadrature (pulls). In spite of the fact that unidentified effects may be still 
affecting some of the data, it is nevertheless interesting to note that, as far as Deuterium 
is concerned, the standard BBN scenario with three neutrinos and WMAP result for 
is quite consistent. 
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Table 2: The main contributions to (cr^) 2 , in percentage, for N e ff = 3.01 and Ub = 0.023. 
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Figure 6: Left plot: the 68 and 95% CL. contours in the cjj, — N e ff plane for WMAP 
(solid lines). The dashed lines are two arbitrary levels of the BBN Deuterium likelihood. 
The joint 95% CL. CMB + BBN region is shown as the filled area. Right plot: the joint 
95% CL. contour of the left plot (solid line) is compared to the same contour with a total 
error on Xd improved by a factor 2 (filled area). 

To analyze this point more quantitatively we have performed a likelihood study using 
the WMAP data and BBN Deuterium abundance only. The likelihood function for D 
has been constructed as described in Section 2. The N e ff and Ub dependence of cr^ has 
been fully taken into account. The result is shown in Fig. |H1 (left plot), from which we 
see the good agreement of a standard BBN with three active neutrinos with the WMAP 
determination of u;&. The product likelihood CMB+BBN is peaked at the values Ub = 
0.023 and N e ff = 3.4. Likelihood marg inalization gives, at 95% CL., u b = 0.023±g;gg| 
and N eff = 3.6l|l, which give at la X D = (2.74i£ff) • 10~ 5 , Y p = 0.256±g;8if and 
X 7Li = (4.7±i;|) • iO -10 - As expected, N eff is not severely constrained either by Xd or 
from CMB: with present uncertainties on Deuterium, the 95% CL. interval corresponds to 
1.3 < N e ff < 6.1. Previous values, when comparisons are possible, are in nice agreement 
with the results reported in [2UI2E]- A more stringent bound on N e ff could be provided 
by A He mass fraction. We will discuss this in more details in the next Section. It is 
however interesting to see whether the value of N e ff could be more precisely fixed by 
Deuterium+WMAP if new and more precise measurements from QSO were available. We 
show in Fig. |H1 (right plot) the results of a simulation with an error reduced by a factor 
two, \J (#22) 2 + (<J2 Xp ) 2 — 0.2 • 10~ 5 . In this case, repeating the same analysis, we obtain 
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2.3 < N e ff < 5.7, still a broad range, but not much wider than what is obtained using a 
customary BBN likelihood analysis of both D and 4 He. 

5 BBN(Deuterium+ 4 J7e) analysis 

The theoretical estimates for 4 ife mass fraction Y p and 7 Li number density Xju show 
a slight disagreement with the corresponding experimental results. For N e ff = 3.01 and 
using the WMAP result for u b we get Y p = 0.2483tg;gg§ and X 7Li = (4.9±H) ■ 10~ 10 . The 
source of uncertainties for these two nuclei is as follows. For Y p the propagated uncertainty 
due to nuclear rates, 044 = 0.0003, is of the same order of magnitude of the one due to 
the variation of uj^. Notice how the smaller error on neutron lifetime, Ar n = 0.8 s, has 
strongly reduced the effects on (cr^) 2 of weak n «-> p processes. Concerning the X 7Li 
uncertainty, this is mainly due to the 4 if e + 3 if e <-> 7 + 7 Be rate, with out, contributing 
for ~ 50% to the total error on X^u- We show in Table 01 the relative contribution of the 
most important rates to (cr^) 2 and (cr^) 2 . 

The status of the primordial abundance measurements for 4 if e and 7 Li is still quite 
involved. There are different determinations of Y p , Y p = 0.234 ± 0.003 jSOj (low Y p ), 
Y p = 0.244 ± 0.002 [SH (high Y p ), and Y p = 0.2421 ± 0.0021 [52J. The latter is the 
result of a recent analysis where the authors performed an accurate study of systematics. 
If the statistical error were not underestimated, these values would be only marginally 
compatible. Waiting for new experimental data, we adopt in our analysis a conservative 
point of view by using 

Y p = 0.239 ± 0.008 . (5.1) 



rate Rk 




rate Rk 




71 <r-¥ p 


29 


4 He + 3 He «-> 7 + 7 5e 


68 


D + D <-> n + 3 He 


29 


3 He + D ^ p + 4 iie 


9 


D + D <-> p + 3 H 


26 


D + £> <-> n + 3 He 


8.5 


p + n <-> 7 + D 


4 


7 Be + n^ 4 He + 4 iie 


7 



Table 3: The main contributions to ((T44) 2 and (cr^) 2 , in percentage, for N e ff = 3.01 and 
Ub = 0.023. n <-> p denotes the weak processes. 
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low 4He 



high 4He 
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-5 -4 -3 -2 -1 

Figure 7: The pulls (see text) of Y p measurements with respect to the theoretical prediction 
for N eff = 3.01 and u b = 0.023, in units of ((crf 4 ) 2 + (af p ) 2 ) 1/2 . 

As it is evident from Fig. where we report the pulls with respects to expected 

1 /2 

theoretical Y p in unit of the combined error ((cr^) 2 + (al xp ) 2 ^j , even the higher value 
of Y p [HI] appears in moderate disagreement (2a effect) with standard BBN. If we let 
N e ff to be a free parameter it is therefore not a surprise that Y p favors a slightly lower 
value for N e ff. In fact, a way to decrease Y p is to delay the freezing out of weak rates 
converting neutron and protons, and this can be achieved by reducing the number of 
relativistic species contributing to the Hubble parameter. In our opinion much effort 
should be devoted to have a more clear understanding of the experimental results or to 
undertake a new measurement campaign for 4 He. On the theoretical side, in fact, we may 
say that the value of Y p is quite robust, and affected by a very small total error. New data 
would eventually tell us whether the standard BBN is fully consistent, or more exotic 
scenarios should be invoked. We discuss one of these, a degenerate neutrino distribution, 
in next Section. 

Concerning 7 Li experimental measurements jSS]-jSZI, the primordial origin of the Spite 
plateau has been recently questioned. In particular the authors of [HI] found evidence for 
a dependence of X-m on metallicity. In any case, there are some indication for a depletion 
mechanism of 7 Li observed in POP II metal poor stars. However, as can be seen in Fig. [HI 
the discrepancy between the most recent observations of |S3]-jnZj and our theoretical value 



17 



7 Li 

Bonifacio et al 2002 
NGC6397 

Bonifacio 2002 M92 

Ryan, Norris & Beers 
'99 

Bonifacio & Molaro 
'97 

-4 -3 -2 -1 

Figure 8: The pulls (see text) of X^u measurements with respect to the theoretical 
prediction for N eff = 3.01 and u b = 0.023, in units of ((af 7 ) 2 + (af p ) 2 ) 1/2 . 

is now reduced to a less than 3a effect. The average over the different observed values for 
X 1L i H2H-E3, which are mutually compatible, gives X 7Li = (1.70 ± 0.17) x 10" 10 . 

As in previous studies we do not use X^u in our BBN likelihood analysis, waiting for 
a more clear understanding of a possible depletion mechanism. Moreover, further efforts 
should be also devoted in trying to reduce the large theoretical uncertainty which is still 
affecting X 7 ii, mainly by lowering the uncertainty on the leading reaction 4 He + 3 He <-> 
7 + 7 Be, presently of the order of 18%. As we mentioned, this produces a large fraction 
of the overall 20% error on Xjn. 

Using the Xq result of [I] and Eq. (|5.1J) for Y p , we consider the likelihood function as 
defined in ()2.4|) . as function of N e ff and Ub- Our results for 68 and 95% C.L. contours 
are shown in Fig. |2l where for comparison we also report the 95% C.L. contour of the 
CMB + D analysis. The effect of 4 ife has been to shift both N e ff and Ub towards 
smaller values. In particular, after marginalization, we get at 2a, uJb = 0.02lto;QQ4 and 
N e jf = 2.5^9. Notice how, with present determination of Y p and X D , the range obtained 
for Ub is rather broad compared with WMAP result. 

We close this section with the theoretical estimate on 3 He and 6 Li number fractions. 
For uj b = 0.021 and N eff = 2.5 we get X 3He = (1.02igi?) • 10" 5 , which saturates the 
recent estimate reported in jHHj for the upper limit of the 3 He to H number density 
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Figure 9: The 68 and 95% CL. contours for the D + 4 He likelihood function in the 
u h — N e ff plane. We also show the result of the CMB + D analysis of Fig. El (left plot). 

from observations of H II regions and planetary nebulae, X^He = (1-1 ± 0.2) ■ 10~ 5 . The 
6 Li primordial density is predicted to be quite small, Xqh = (1.2 ± 1.7) • 10~ 14 , which 
results in a ratio 6 Li/ 7 Li = (2.8 ± 4.0) • 10~ 5 , much smaller than the experimental value 
6 Lip Li = 0.05, measured in two stars in the galactic halo (see, for example, [59J and 
refs. therein). Very close result are obtained for the other analysis considered in the 
previous Section. We recall however that the theoretical estimate is still suffering the 
large uncertainty affecting the rate A He + D «-> 7 + 6 Li, which is supposed to give the 
main production channel during BBN (see our discussion in Section 2). Actually any 6 Li 
detection has important consequences in constraining the destruction rate of 7 Li in POP 
II stars, as well as exotic processes to D synthesis, as for example photodissociation or 
spallation of 4 He jHO] • A summary of the ranges for Ub and N e f f and of nuclei yields is 
reported in Table 4. 

6 Bounds on neutrino chemical potentials. Still room 
for extra relativistic species 

The cosmological effects of non zero neutrino antineutrino asymmetries have received 
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CMB + BBN(D) 
at N eff = 3.01 


CMB + BBN(D) 


BBN(D + 4 He) 


u b (2 a) 


n noo+0.003 


n23 +0 - 003 


n a 91 +0.005 


N eff (2 a) 


3.01 


3.6±1;| 


O r+1 1 
2-5-0.9 


X D -10 5 (la) 
(2.78±°;||) 


o c:c:+0.22 
Zj.oo n la 


2.74 + , f 

— U.Ol 


9 7 c-+0.53 

— U.Ol 


Y p (la) 
(0.239 ± 0.008) 


0.2483+O™ 


— U.U10 


0.240 ±0.008 


X 7Li .10 10 (l<r) 
(1.70 ± 0.17) 


4.9±\ A 2 


4.7±i:i 


4.3±i : g 


Xz He -ltf(la) 
(1.1 ±0.2) 


n QQ+ - 07 

u - yy -0.08 


1 m+ - 09 

i - u± - 0.10 


1 02 +0 - 10 
i - uz -0.11 


XeLi/X 7Li (la) 
(0.05) 


(2.4 ±3.4) • 10~ 5 


(2.8 ±3.7) • 10~ 5 


(2.8 ±4.0) • 10~ 5 



Table 4: A summary of the results obtained in the three different SBBN likelihood analysis 
discussed in this work. For nuclear abundances we report in parenthesis the experimental 
value or the best estimate currently available: actually, only in some cases a direct com- 
parison is possible, as explained in the text. The larger uncertainties affecting Y p in the 
third and fourth columns are due to the variation of N e ff in the 1 — a range. 
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a renewed attention in the last few years (23] > The bounds from BBN and CMB 

data on the asymmetry parameters £ a for the three active neutrinos, a = e, fi, r, are more 
stringent by considering the effect of oscillations in primordial universe. In fact, in [33] it 
was shown that for the large mixing angle solution to solar neutrino problem, presently 
favored by data, oscillations tend to homogenize the which are therefore all bounded 
by the more stringent constraint on £ e . Neutrino chemical potentials affect BBN via their 
contribution to energy density, 



and so to Hubble parameter. In particular <^ jT only enters in changing the value of N e fj. 
On the other hand £ e also affects the n/p chemical equilibrium at freeze out, changing the 
neutron to proton ratio as exp(— £ e ). A positive £ e (more neutrinos than antineutrinos) 
favors n —>■ p reactions with respect to inverse processes, and this eventually translates 
into a smaller value for Y p . 

In this section, using the WMAP result on Ub as a prior, we use our BBN code to 
update our previous analysis |2| of degenerate BBN (DBBN), with two main purposes: 

i) since, as we stressed in the previous Section, the Y p experimental value is slightly 
lower than what expected in standard BBN, we analyze how the agreement im- 
proves allowing for non zero £ a , but fixing the number of relativistic species, besides 
photons, to the three active light neutrinos. 

ii) in the general case we determine the upper bounds in DBBN on N e ff and £ a . 

We construct the BBN likelihood function C BBN (N e ff,uJb,^e) defined in Eq. ([2.4)1 . 
with the extra dependence on the parameter £ e . The use of Y p may seem inconsistent 
with our discussion of the previous Section. For standard BBN in fact we stressed that 
the experimental value of Y p shows a slight disagreement with the theoretical estimate, 
possibly a sign of unaccounted systematic effects. However the introduction of a (positive) 
£ e in the DBBN case leads to a smaller theoretical value for Y p even for N e ff ~ 3, so in 
this case one may interpret the experimental findings as an indication of a neutrino- 
antineutrino asymmetry. 
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Figure 10: The 68 and 95% contours for the D + 4 He degenerate likelihood function in 
the uj b — £ e plane in the scenario with only three active neutrinos and without prior on ou b 
(left plot). Right plot shows the marginalized likelihood versus £ e when a prior on uo b is 
imposed, u b = 0.023 ± 0.002. 

In the minimal scenario we assume only three active neutrinos contributing to N e ff. 
Since £ MiT are enforced by oscillations to be of the same order of magnitude of £ e , we can 
write down a constraint for N e ff, 



N, 



eff 



3.01 + 3 



30 


(I)' 


15 




y 




- 





NeffiQ ■ 



(6.2) 



In the second scenario we parameterize possible extra light degrees of freedom with AN e f f , 
defined by 

Neff = #.//(&) + AiV e// . (6.3) 

In both these cases we adopt the experimental values/errors for Y p and Xr> reported in 
Section 4 and 5. 

We show in Fig. ^] (left plot) the 68 and 95% CL. contours for the first scenario. When 
imposing the WMAP prior with a conservative error (uj b = 0.023 ± 0.002) we obtain the 
marginalized likelihood for £ e reported in the right plot of Fig. giving £ e = 0.03lo;ocs 
(2a). Correspondingly X D = 2.51 • 10~ 5 , Y p = 0.2414. 

Fig. 1111 summarizes our findings in the second scenario. At 95% C.L., we get N e ff = 
3.9^314 and £ e = 0.131q;26, again using the WMAP prior on u b . We notice that our limits 
on N e f f are broader than what found in , since we use only the uj b prior from CMB in 
the analysis. This results in a more conservative bound, independent on the assumption 
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Figure 11: The 68 and 95% CL. contours in the £ e — AN e ff plane for the degenerate 
BBN likelihood function. 

that N e ff is the same at both the CMB and BBN epochs. Our upper limit on N e ff, 
being sensibly larger than 3, should be understood as a direct bound to possible extra 
light particles. Notice that in the DBBN scenario a fourth sterile neutrino, as for example 
required to interpret LSND evidence for <-> V e oscillation, is not yet ruled out. See also 
for a neutrino mass analysis of this issue. 

7 Conclusions 

In this paper we have discussed the present status of Primordial Nucleosynthesis after 
the recent measurements of CMB anisotropies by the WMAP experiment, which repre- 
sents a new relevant step towards a better understanding of the evolution of the universe, 
and in particular gives a very precise determination of the baryon density parameter Uf,. 
The question of whether the standard picture of BBN is a satisfactory scenario for produc- 
tion of light nuclei requires an increased precision in both its theoretical and experimental 
aspects. 

As discussed in the paper, it seems to us that some results in this direction have already 
been achieved. The accuracy of theoretical estimates has been improved by refining the 
estimate of n/p ratio at decoupling. Most importantly, a complete reanalysis of all nuclear 
rates which enter the BBN network has been made, using the NACRE compilation, as 
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well as all most recent available results. This analysis has been used to implement a new 
BBN numerical code. Particular attention has been devoted to the rates which mainly 
contribute to evolution of observed nuclei abundances. Moreover several processes which 
were not considered in previous analysis have been added. The main outcome of this study 
is a reduced theoretical uncertainty of the BBN model, which is presently well below the 
corresponding experimental errors, at least for D and A He. The 7 Li number fraction is 
still affected by a 20% uncertainty, mainly due to the effect of A He + 3 He <-> 7 + 7 Be rate. 
More accurate measurements of this rate are therefore strongly desirable. 

From the experimental point of view, Deuterium determinations show a very good 
agreement with the corresponding theoretical expectation, and are fully compatible with 
the value of u>j, suggested by WMAP data. This is quite remarkable, since D is strongly 
dependent on the baryon content of the universe. On the other hand we do think that new 
measurements of 4 if e and 7 Li should be performed, along with a clear understanding of 
possible depletion mechanisms of primordial X-jh. It should be mentioned however that 
the most recent results for 7 Li, compared with our theoretical prediction, reduce the 
disagreement for X^n at less than 3a level. The A He mass fraction still has two distinct 
determinations which are mutually incompatible. Even adopting a conservative value, as 
we did in this paper, still the amount of A He experimentally detected seems slightly lower 
than expected for a standard BBN scenario. This may be seen either as due to possible 
systematics in the measurements and/or extrapolation technique for Y p , or rather as the 
effect of exotic physics in the BBN scenario, as for example neutrino degeneracy. For an 
alternative method, still affected by large uncertainty, to infer Y p see Ref. jHEl- Possible 
measurements of other nuclei, like 3 He, for which there is a single recent estimate, or 
6 Li, though difficult 1 , may greatly help in clarifying the overall soundness of Primordial 
Nucleosynthesis scenarios. 

The main results of this paper are reported in Table 4 while our analysis can be 
summarized as follows: 

i) we have described some key examples of the study performed to improve the accu- 

1 Actually, only for a smaller uib, a large 4 He + D — > 7 + 6 Li reaction rate, at the upper limit of the 
present uncertainty range, and almost no depletion mechanism in act in hot and metal-poor PopII halo 
stars, one could reveal the primordial (as opposed to the cosmic-ray produced) 6 Li in the near future. 
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racy of the nuclear network adopted in our code. A full description will be presented 
elsewhere; 

ii) we have performed a likelihood analysis of the WMAP results and BBN (Deuterium) 
abundance on u b and N e ff, showing that the simplest standard BBN, with three 
neutrino as the only relativistic particles besides photons at the BBN epoch is largely 
consistent. In particular we found, at 95% C.L. uo b = 0.023±o;qo2 an d N e ff = 3.6tl'^- 
Notice the weak sensitivity of both D and CMB data to N e ff. We stressed however 
that an improved accuracy of D measurement may provide more accurate bounds 
on this parameter, independently of information on Y p ; 

hi) we have studied, again via a likelihood method of the Xp and Y p observables, 
the theory versus experiment status. We find that taking into account the present 
determination of A He shifts the preferred values of both N e ff and uj b towards smaller 
values than in the CMB + D analysis, u b = 0.021^;^ and N eff = 2.5±^ at 95% 
C.L.. Finally we have also reported our estimate for 3 He fractional density and 
6 Lip Li ratio (see Table 4); 

iv) in the framework of degenerate BBN we reported the bounds on possible extra rel- 
ativistic degrees of freedom. Oscillations enforce the three active neutrinos to share 
the same small chemical potential at the BBN epoch and we find as a common 
bound —0.03 < £ e ,/i,T < 0.11 at 95% C.L. for N e ff corresponding to three active 
neutrinos. In this case the agreement of A He theoretical expectation with experi- 
mental data improves. If we leave N e ff as a free parameter we get N e jf = 3.9^4 at 
95% C.L.. In particular the upper bound represents the largest extra contribution 
to the relativistic energy density due to exotic particles consistent with BBN. 

After the huge progresses in the study and detection of Cosmic Microwave Background 
anisotropies, we think that further efforts in refining the BBN theoretical and experimental 
aspects may provide new pieces of information on cosmology and fundamental physics. 
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